Lossy Quantum Optical Metrology with Squeezed States 
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We study the precise phase estimation using squeezed states with photon losses present. Our exact 
quantum Fisher information calculation shows significant quantum enhancement and thus reveals 
the benchmark for practical quantum metrology in this noisy scenario. However, we find that the 
existing parity measurement scheme [P.M. Anisimov et al, Phys. Rev. Lett. 104, 103602 (2010)] 
behaves worse than even classical cases given very small losses, unless we take an appropriate loss 
dependent phase shift. Using our formulae, the two optimized aspects including the pre-detection 
phase shift and the consequent light intensity of a tradeoff strategy for photon resource arrangement 
can be both calculated. Therefore our result makes it possible to experimentally realize quantum 
metrology of phase estimation with squeezed states. 



INTRODUCTION 

Many important tasks in scientific research involve 
physical processes that contribute to phase measurement. 
Precice measurement is of nontrivial significance since 
physical laws (e.g. the uncertainty principle) impose lim- 
its on it. By finding measurements that optimally resolve 
neighboring quantum states, Braunstein and Caves [l| 
proposed a metric on density operator space and a more 
general uncertainty principle, which lead to a bound 
S(/> > 1/ \/vFq of the minimal error while estimating a 
parameter (/>, where Fq is the quantum Fisher informa- 
tion(QFI) of the probe state used in an experiment and 
v is the number of identical experiments carried. Based 
on [l[ , by employing quantum entanglement Q , quantum 
metrology promises the possibility of attaining a more 
fundamental Heisenberg limit (HL) 1/n (n is the number 
of particles involved) rather than the standard quantum 
limit(SQL) 1/y/n for precison of conventional phase esti- 
mation. Thus the knowledge of QFI has become the key, 
and by comparing QFIs of different states we can specify 
the probe state to beat the SQL. Measurement strat- 
egy is another challenge for researchers, since within the 
framework we can see that quantum metrology requires 
suitable choice of measurement basis. Many schemes up 
to now originate from the mathematically derived POVM 
measurement for pure state [1] whose physical realization 
remains unclear. 

While the framework is appealing, imperfections in 
practical realizations result in certain amount of dete- 
rioration of the quantum enhancement . Besides, an 
analytical form of the exact QFI is formidable to ob- 
tain due to the difficulty in diagonalizing the compli- 
cated density matrices after noisy evolution. Representa- 
tively, specific [H or general [Hj upper bounds and numeri- 
cal analysis [3| of QFI with noises have been put forward. 
What's more, the real- world noise makes the measure- 
ment task more complicated and intractable as well. For 
instance, we find in this paper that a good scheme for 



phase measurement of squeezed vacuum state will be- 
have quite worse than even classical cases once photon 
losses appears. 



Aiming for a deeper understanding of this malfunction 
of measurement scheme, we study noisy phase estima- 
tion with mathematical rigor for quite general squeezed 
Gaussian states in this paper. We firstly analyze the en- 
tangled two-mode case with exact QFI for lossy channels 
calculated via the fidelity approach, and the reason for 
the singling out of two mode squeezed vacuum(TMSV) 
state is also given. Therefore, much more precise mastery 
of squeezing for practical quantum optical metrology is 
achieved, which still show significant quantum enhance- 
ment. In order to exploit this enhancement, we secondly 
carefully analyze the behavior of parity measurement, 
which was originally proposed for trapped ions[8| and also 
appeared in pi] [9] , applied on our aforementioned prefer- 
able TMSV state with photon losses taken into account. 
While giving explicitly the exact phase measurement er- 
ror to show the malfunction, we address on our newly 
discovered optimal phase measurement point depending 
on noise, ie. knowledge of channel losses helps to deter- 
mine where to measure. This leads to vital modification 
of ordinary measurement procedures for ideal cases, and 
implies some general properties for noisy cases. For in- 
stance, unlike the ideal scenario, a balance between the 
average photon number n used in a single experiment 
and the repeated time v of identical experiments should 
be built if total photon resource provided. The contem- 
porarily realizable measurement scheme we devised to 
adapt for these new features arisen in noisy metrology 
turns out to have comparable improvement than classi- 
cal setups. Thus, our exact QFI evaluation and measure- 
ment scheme as a whole undoubtedly justify the quantum 
enhancement in practical optical metrology, which could 
be heuristic and instructive for theorists and especially 
for experimentalists. 
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FIG. 1. Setup for quantum metrology with MZI 



EXACT QFI WITH PHOTON LOSSES 



The schematic display of the two-arm interferomet- 
ric quantum metrology is shown in FigJT] The essential 
method of phase estimation in quantum optical inter- 
ferometry consists of light's passing through a Mach- 
Zehnder interferometer (MZI) and measurement at the 
output ports. A typical MZI consists of two 50:50 beam 
splitters and a phase shifter in one arm. The descrip- 
tion of noise involves introducing ancillary modes of the 
environment. Fictitious beam splitters with transmis- 
sivity 771,772 are used to produce photon losses. The 
losses accompany with phase accumulation simultane- 
ously. But this can be equivalently simulated by commu- 
tative seperate phase shifting and photon losing, which 
can be easily seen from the fictitious beam splitter's 
Kraus operator form for the dispersive arm's mode. For 
a Gaussian state [lOj, the initial covariance matrix is 
7o © Ie with possible displacement. The state after pho- 
ton losses and phase accumulation is characterized by 

7 is 



M^M loss ( l0 ® IeW^M* ={jr where 

Mioss/^ represent the linear transformation of field op- 
erators due to the devices. Leaving out the environ- 
ment, the covariance matrix for the system we consid- 
ered is 7 = 7(0). Firstly, contrary to displacement op- 
eration's uselessness for quantum metrology due to its 
pure translation of variance ellipse in quadrature space, 
squeezing notably changes the variance within the limit 
of uncertainty principle of the two conjugate quadratures, 
which strongly implies its advantage in quantum metrol- 
ogy. In addition, simple physical intuition implies that 
finite temperature degrades the quantum coherence of 
a state. Lastly, compared with one-mode counterparts, 
two-mode Gaussian states have more complex and in- 
triguing properties, and most importantly, they are com- 
pulsory in MZI parity measurement (see MZI part). Thus, 
we concentrate on TMSV states for lossy interferome- 
ters rather than displaced two-mode squeezed thermal 
(DTMST) states. We point out that the mathematical 
derivation of the general DTMST states notably resem- 
bles the following TMSV case. Besides the change of the 
aforesaid 70 matrix, the two cases' covariance matrices 
after losses share the similar form like (pQ), and this is 
central for our discussion. 



Calculation shows that for an initial TMSV state 
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1 + 2772 sinh r, 
squeeze param- 



where d\ = 1 + 2771 sinh r, d 2 
a = ^r]i 772 sinh 2r and r is the 
eter. Our general QFI calculation method relies 
on the relationship between the Bures distance and 
QFlQ: F Q = 4[ds B ures(p + d p,p)] 2 /d0 2 , where 
ds Bu res(pi,P2) = \/2[l - F(p x ,p 2 )] is defined using fi- 
delity F(p u p 2 ) = TrV^ 



/PiP2vPi- Different expres- 
sions for fidelity have been proposed [ill 12]- We uti- 
lize diagonalizable quadratic Hamiltonian transformation 
method [12| for our use. We conclude from (pQ) that the fi- 
nal state after evolution is a two-mode squeezed thermal 
state, whose density operator surely has a diagonalizable 
quadratic Hamiltonian form. Thus, by feeding it to [l2|'s 
method, we get the exact lossy QFI: 



F c 



2n(n + 2)771772 
2 + 71(771 +772 - 2771772) 



(2) 



where n = 2 sinh 2 r is the average photon number of the 
initial ideal TMSV state. Reduction to ideal case Fq = 
n(n + 2) (we call it "modified HL") is straightforward. 
The power of our calculation method manifests itself for 
this generality of dealing with different lossy scenarios. 
For equal losses in both arms (771 = 772 = 77) and one- 
arm losses (771 =77,772 = 1), the QFIs are ^^2^ an d 

n(™-^>+2 • ^ ne l a tt er is more practical since losses are 
mostly generated by the phase accumulation. 

Here we can see our exact lossy QFI (|2|) (corresponds 
to "quantum limit" 1 / ^/FQ in FIG|2j) for both particular 
cases ascertain the general framework in [5[ , where their 
upper bound of QFI for optical interferometry gives scal- 
ing 1/y/n for n ^> 77/(1 — 77) and 1/n for n <C 77/(1 — 77). 
Despite of the same physical model, our intuitive upper 
bound of QFl[l3| generally acts worse than our exact 
QFIs as depicted. This strongly endorses the importance 
of our exact QFI study. As for the performance of TMSV 
state, the QFIs for two lossy scenarios approach and are 
a little worse than their classical limits in very large loss 
limit (i.e. 77 0). Notably, TMSV state beats SQL even 
in the presence of large losses (77 = 0.6 for two-arm losses, 
77 = 0.4 for one-arm losses). This implies that cleverly 
designed measurement scheme is possible to achieve en- 
hancement with squeezing in a large extent of practical 
setups. 



MEASUREMENT SCHEME 

We denote the evolution and the expectation value 
of parity measurement at the output mode respectively 
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FIG. 2. (LEFT) Estimation precision 8(j) vs average photon 
number n. Black dashed- dotted line: Modified HL. Black 
dashed line: SQL. Green and blue line: Quantum limits for 
two-arm equal- loss (771 = 772 = 0.8) and one- arm- loss (771 = 
0.8,772 = 0) models respectively. Red line: 
(RIGHT) Estimation precision 5<j) vs transmissivity 77. Black 
dashed-dotted line: Modified HL. Black dashed line: SQL. 
Green and blue solid line: Quantum limits for two-arm equal- 
loss and one-arm-loss models respectively. Green and blue 
dashed line: Classical limits^ l/y / nrj for two-arm equal-loss 
and for one-arm- loss models respectively. Red line: [l3[. 



as |^/) = Ui OSS yMZi \^TMSv) ah ® |0,0) e and (II) = 

(^/|n|^/), where UiossyMZi = UbsU^Uioss&bs- Us- 
ing exponential operator reordering method [14] hencefor- 
ward, we can calculate both the UiossyMZi operator and 
expectation values of observables. Note that they are 
not trivial to obtain. For instance, it's easy to show that 
two one mode squeezed vacuum product state with same 
average photon number in each mode passing lossy MZI 
carries no phase information within parity measurement. 

Parity operator II = exp (— ma) a) represents the ob- 
servable to be measured. Research shows that parity 
measurements enables sub- SQL interferometric phase es- 
timation in a wide range of quantum photonic states [9] 
with no need to know the complete photon counting 
statistics. Notably, it's useful for both qubit-like very 
nonclassical states (e.g. NOON state) and continuous vari- 
able quantum states (e.g. our indefinite photon number 
TMSV state). Besides, its unambiguous physical mean- 
ing and realization stand out among other mathematical 
measurement basis. It's in this spirit we devote to exploit 
the quantum enhancement via parity measurement. 

Some concern about parity measurement is due to the 
absence of highly reliable photon number resolving detec- 
tor, which seems to be necessary for ITs dependence on 
photon number operator. However, parity measurement 
actually could be done without photon number resolv- 
ing detectors. We highlight the equivalence between the 
value of Wigner function at the origin and the parity of 
one mode field [lit: (ft) = ^W(0,0). Quantum tomogra- 
phy techniques provides practical construction method of 
the Wigner function of a radiation field. The fact that we 
need W(0, 0) alone leads to an explicit experiment setting 



FIG. 3. (LEFT) Optimal measurement point cj) vs average 
photon number n and transmissivity 77. 

(RIGHT) Phase estimation error 5(j) vs average photon number 
n for one-arm losses (77 = 99%). Purple line: Classical limit [jj. 
Blue line: Quantum limit of exact QFI. Green line: Optimal 
phase measurement error via parity detection. 



using balanced homodyne detection [16j, wherein only a 
few extra common optical devices and standard intensity 
difference measurements are necessary. We consider this 
method based on highly developed detector technologies 
very favorably for a demonstration of our scheme. Very 
clear and detailed discussion of the setups and mathe- 
matical derivations can be found in [16j. 

Our calculation gives the measurement result and the 
variance of phase estimation that quantifies the quantum 
enhancement explicitly respectively as 



= 



<n>« 



<n 2 >, - (n) d 



2' 



w(nw + 2) s 



(d(ii) 



2rj{ t]\n(n + 2) 2 sin 2 2f 



(3) 



(4) 



where uj = 2f]if] 2 (n + 2) cos 2 (j) + (771 + 772) (2 - 771 - ^2)- 
dUdl]) reduce to the ideal case [7] when 771 =772 = T 
(ft) , = 1 = i+"("+2)cqb 2 ^ Thig ideal 

N l( P ^1+^(^+2) cos 2 r I sin 01^/71(71+2) 

Acj) is a monotonic function of <j) if <j) G [0,7r/2], i.e. the 
best precision is achieved in the vicinity of <fi = tt/2. How- 
ever, we can notice fl3J)'s explicit dependence on both n 
and 77, i.e. the monotonic property against (j) no longer 
exists in the lossy case. Actually, this new relationship 
A(j) = A0(n, 771, 772, (j)) is not simple. 1) It will be strongly 
affected by the values of n and 77, where the average pho- 
ton number n directly depends on the squeeze parameter 
£ of TMSV state and 771 , 772 is due to the degree of im- 
perfection of the phase shifter or the arms. 2) Instead of 
exploitation of quantum enhancement at the vicinity of 
(j) = 7r/2, we can find an optimal phase <p to measure. 
This optimal phase <p = (j) (n, ?7i, ^2), whereupon we at- 
tain the minimal A0, can be straightforwardly calculated 
by simple derivative from our exact (|4]). 3) We point out 



that without this "optimal measurement point" amend- 
ment, the direct phase estimation will perform very badly 
and no sign of enhancement could be observed (so we 
don't show it in FIG BRIGHT). This is indeed a quite 
unexpected feature in the presence of noise, and no trace 
can be seen in the ideal case. Since losses are mostly 
generated by the phase accumulation, we focus on the 
Vi = V2 =0 case henceforth. Interestingly, there do ex- 
ist monotonic dependences on n and r] respectively of this 
O , which is shown in FIGlSfLEFT). Thus, we found that 
the salient and quite notable difference between ideal case 
and practical lossy case is the varying optimal measure- 
ment point, which will obviously affects the real- world 
experimental strategies for quantum metrology. Physical 
considerations, e.g. the generality of quantum metrol- 
ogy framework and the mathematical equivalence be- 
tween quantum metrology in optical system and some 
other physical systems, imply that the above discussion 
won't be limited to specific state preparation, phase ac- 
cumulation or detection scheme and these new features 
of noisy quantum metrology will appear extensively in 
many other cases. 

Compared with our exact lossy QFI and the classi- 
cal limit yjt= for coherent state [3|, the performance of 
parity measurement at optimal measurement point with 
losses in the dispersive arm alone has been calculated, 
where significant quantum enhancement is achieved in 
r] = 0.99 case as shown in FIG RIGHT). To our knowl- 
edge of contemporary optical technology in laboratories, 
this is already enough for a proof-of-principle demonstra- 
tion of our scheme . In addition, our exact lossy QFI 
bound of phase estimation cannot be totally achieved 
with parity measurement since QFI provides the tight- 
est bound of mixed state after photon losses within the 
framework of quantum metrology, i.e. there's still more 
effort to be devoted for the quantum extraction of phase 
information. Although being monotonic, our parity mea- 
surement precision slowly deviates from the lossy QFI 
with large average photon number and imperfection, this 
proves the vulnerability of quantum enhancement against 
noises. 

In practical experiment settings, we have to make bi- 
ased measurements in order to achieve the best precision. 
As aforesaid, we firstly calculate the optimal phase <p 
from the knowledge of state preparation and experimen- 
tal imperfection. If the phase (j) to be measured is small, 
we then adjust the phase shifter to a biased phase shift 
(f = cj) + (j). The parity measurement of cp shares a same 
precision with 0, which can be easily seen from (j4]). In 
real need of precise phase measurement, we can always 
make the phase to be measured small through the knowl- 
edge of specific phase accumulation physical processes. 

As for the situation closer to real experimental setups, 
we generally perform the experiment of identical setting 
many times, i.e. we divide our N-photon resource into 




FIG. 4. Repeated-experiment phase error 5<j) vs average pho- 
ton number n for one- arm losses (77 = 99%, 98%, 97%, 96% 
from left to right, N=200 fixed). Red horizontal line: Clas- 
sical limit [3]. Blue line: Quantum limit of exact QFI. Green 
line: Optimal phase measurement error via parity detection. 



v identical quantum states of light with n = N/ v being 
interpreted as the average photon number in the above 
discussion. Therefore the more general formula of mea- 
surement error is 5<fi(n : 77, = ^^==^-. With N fixed, 

the optimization of phase measurement point is identical 
to (j4j). However, there turns out to be an optimal n Q 
that increases with 77 to attain minimal error Scj) for fixed 
77 and N as shown in FIGHJ i.e. very large photon num- 
ber is not simply always good, which means there exists 
a tradeoff between the quantum enhancement surpassing 
the classical error reduction and the more losses 

inevitably accompanying the more photons. Therefore 
noisy scenario requires a new arrangement of total pho- 
ton resources, since the quantum enhancement in fact has 
a very different asymptotic performance than the simple 
\j \fv or ideal HL 1/n. 
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